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2^ ■ Abstract 



The WDW equation of arbitrary Hartle-Hawking factor ordering for several minisuperspace universe 
models, such as the pure gravity FRW and Taub ones, is mapped onto the dynamics of corresponding 
classical oscillators. The latter ones are studied by the classical Ermakov invariant method, which is 



> 

o 
o 

' a natural aproach in this context. For the more realistic case of a minimally coupled massive scalar 

field, one can study, within the same type of approach, the corresponding squeezing features as a 

CN ■ 

^ . possible means of describing cosmological evolution. Finally, we comment on the analogy with the 

00 ■ accelerator physics. 

O 

OO ' 1. The formalism of Ermakov-type invariants [1] can be a useful, alternative method of 

investigating evolutionary and chaotic dynamical problems in the "quantum" cosmological 
! framework [2] . Moreover, the method of adiabatic invariants is intimately related to geomet- 
Q \ rical angles and phases [3] so that one may think of cosmological Hannay's angles as well as 
^\ various types of topological phases as those of Berry and Pancharatnam [4]. In the following 
we apply the formal Ermakov scheme to some of the simplest cosmological pure gravity oscil- 
lators, such as the empty Friedmann-Robertson- Walker (EFRW) "quantum" universes and 
the anisotropic Taub ones. The EFRW Wheeler-DeWitt (WDW) minisuperspace equation 
rN \ reads [5] 



where Q, assumed nonzero [6] and kept as a free parameter, is the Hartle-Hawking (HH) 
parameter for the factor ordering [7], the variable Vl is Misner's time [8], and n is the 
curvature index of the FRW universe; k = 1,0, —1 for closed, flat, open universes, respec- 
tively. For K = ±1 the general solution is expressed in terms of Bessel functions, = 
(Ci/,(ie"2^) + C2i^a(ie-2^)) and = e'^"^ (Ci J,(ie-2^) + CaF^de-^^)), re- 

spectively, where a = Q/4:. The case k = is special/degenerate, leads to hyperbolic 
functions and will not be dealt with here. Eq. (1) can be mapped in a known way to the 
canonical equations for a classical point particle of mass M = e'^^, generalized coordinate 
g = \Ef, momentum p = e'^^\E', (i.e., velocity v = and identifying Misner's time Q with 
the classical Hamiltonian time. Thus, one is led to 



(2) 



(3) 



These equations describe the canonical motion for a classical EFRW point universe as derived 
from the time-dependent Hamiltonian of the inverted oscillator type [9] 

i/ei(^^)=e-'3"^-/^e(«-')^^ . (4) 
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For this classical EFRW Hamiltonian the triplet of phase-space functions Ti = ^, T2 = pg, 
and T3 = ^ forms a dynamical Lie algebra (i.e., H = J2nf^n{^)Tn{p,q)) which is closed 
with respect to the Poisson bracket, or more exactly {Ti,T2} = — 2Ti, {72,23} = — STa, 
{Ti^T^} — —T2. Using this algebra Hd reads 

= e-^^Ti - ne^^-^^^'Ts . (5) 



The Ermakov invariant X belongs to the dynamical algebra, i.e., one can write X 

ax 
an 



Y.r ^ri^)Tr, and by means of ^ = —{I,H} one is led to the following equations for the 
functions er{fi) 



+ E 



er 

n 



0, (6) 



where are the structure constants of the Lie algebra that have been already given above. 
Thus, we get 

= _2e-«^e2 

e2 = -«e(«-4)^^ei-e-«% (7) 
^3 = _2«:e(^-^)"62 . 

The solution of this system can be readily obtained by setting ei = giving €2 = —e^^pp 
and €3 = e^'^^f? + where p is the solution of the Milne-Pinney (MP) equation [10], 

p + Qp — K,e~^^p = ^ There is a well-defined prescription going back to Pinney's note 
in 1950 of writing p as a function of the particular solutions of the corresponding parametric 
oscillator problem, i.e., the modified Bessel functions in the EFRW case. In the formulas 
herein, we shall keep the symbol p for this known function. In terms of the function p(il) 
the Ermakov invariant reads [11] 

Iefrw- ^ + — -^[p^^-p^^) +^[—) ■ (8) 

2. Next, we calculate the time- dependent generating function allowing one to pass to new 
canonical variables for which I is chosen to be the new "momentum" S{q, P = I, e{fl)) — 
P dq p{q , X, e{fl)) leading to 
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S{q,X, e{n)) = e^"|-^ +Xarcsin 
2 p 



where we have put to zero the constant of integration. Thus, 9 — ^ — arcsin(^^ 



Moreover, the canonical variables are now q = py/2Tsm 6 and p = cos 6 + e'^^pp sin 6 

The dynamical angle will be A^'^ = Jl^^{^^)dQ' = Jo[^^ - ^J^(^)]dQ', whereas the 

geometrical angle reads Ad^ = ^ /si^o[^(e'^^ pp) — 2e^^ p'^]dil' . Thus, the total change of 

angle is = Jq^ ^—^ — dfi'. On the Misner time axis, going to — oo means going to the origin 
of the universe, whereas f2o = means the present epoch. Using these cosmological limits 
we obtain the interesting result that the total change of angle A9 during the cosmological 
evolution in Q time can be written up to a sign as the Laplace transform of parameter Q of 
the inverse square of the MP function, A^ = — Li/p2(Q). 

3. We now sketch the minisuperspace Taub model for which the WDW equation reads 

where V{(3) = |(e~^^ — 4e~^^). This equation can be separated in the variables Xi = 
—AO, — 8/3 and X2 — — 4Q — 2/3. Thus, one gets the following two independent ID problems 
for which the Ermakov procedure can be repeated along the hues of the EFRW case 

, Q , 1 -U„ = (11) 



and 



dxi 12 dxi \ 4 144 



dx2 3 dx2 \ 9 



+ (u'--e^')^T2^0 . (12) 



The quantity uj/2 is the separation constant. The solutions are \t'Ti = ^Tai = 
ei-Q/^^ :''^Z,,,^{ic'''-''^/ 6) and ^t2 = 'in., = e^^l^>^ Zia^{i2e''^l'^ /?,), respectively, where 
ai = - (Q/12)2 and as = ^JAu'' - (Q/3)2. 

4. A more realistic case is provided by the minimally coupled FRW-massive-scalar-field 
minisuperspace model. The Ermakov approach, which differs from the previous one, will be 
studied in detail elsewhere. The WDW equation reads 

[dl + Qda -dl- «e-^^ + m2e-^^(/)2]*(Q, 0) = , (13) 

and can be written as a two-component Schroedinger equation (see e.g., [12]). This allows 
one to think of cosmological squeezed states based on the Ermakov approach [13,14]. For 
this one makes use of the factorization of the Ermakov invariant X = ^(66^ + |), where 
h = (2;i)-i/2[2 + i{pp - e^^^pg)] and = {2h)-^/^[^ - i{pp - e^-^pg)]. is a fixed 
HH factor ordering parameter. Let us now consider a reference Misner-timc-indcpcndcnt 
oscillator with the Misner frequency Uq corresponds to an arbitrary epoch Qq for which one 
can write the common factorizing operators a = {2tiuo)~^^^[uJoq+ip], = {2huJo)~^^'^[uJoq— 
ip]. The connection between the a and b pairs is given by b{Q) = p{fl)a + z/(f2)a^ and 
b^{n) = p*{n)a^ + z/*(fi)at , where /i(fi) = {AiUo)-^/^p-^ - ie^^^p + uop] and iy{n) = 
{4ujo)-^/'^[p~^ - ie'^'^^p - ujqp] fulfiU the well-known relationship \p{0)\^ - \iy{0)\'^ = 1. The 
corresponding uncertainties are known to be (Ag)^ = -^Ifi — ^l"^, (Ap)^ = ^\fi + ^^P, and 



{Aq){Ap) — + 7y| l/i — z^l showing that in general the Ermakov squeezed states are not 
minimum uncertainty states [14]. 

5. As was first noticed by K.R. Symon [15], the Ermakov invariant is equivalent to the 
Cour ant- Snyder one in accelerator physics, which defines the admittance of the accelerating 
device. This allows in a certain sense a beam physics approach to cosmological evolution. 
The point is that under the assumption of no coupling between the radial and the vertical 
betatron oscillations, the latter ones are described by the Hill equation z + n{s)KgZ = 0, 
where n is the magnetic field index, and Ko is the curvature of the orbit which is parametrized 
by s that may be considered as a counterpart of fl. The solutions can be written as z± ~ 
w(s)e='='^^*^ and only a real '4> leads to bounded oscillations. The amplitude w(s) satisfies 
a MP equation and moreover ip = w^^. However, in "quantum" cosmology, ip is a pure 
imaginary action functional leading to instabilities in the \E' solutions. In other words, while 
in accelerators we are interested in stable periodic solutions, in "quantum" cosmology there 
are the unstable parametric solutions that come into play. 
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